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We  have  used  an  ab  initio  SCF  molecular  orbital  approach  in  conjunction  with  the  isodesmic  reaction  procedure  to  investigate 
anomalous  energy  effects  in  aliphatic  and  aiicyclic  aza  systems.  These  molecules  are  found  to  be  significantly  stabilized, 
which  we  attribute  to  o-conjugation  of  the  nitrogen  lone  pairs;  these  are  presumably  delocalized  to  some  extent  through  the 
o-bond  framework.  The  stabilization  increases  with  the  number  of  nitrogens  in  a  series  of  similar  molecules  and  is  greater 
in  cyclic  systems,  particularly  when  rings  are  fused.  The  magnitudes  of  the  most  negative  electrostatic  potentials  of  the  aza 
nitrogens,  which  can  be  viewed  as  indicating  the  degrees  of  delocalization  of  the  nitrogen  lone  pairs,  decrease  as  the  added 
stability  increases.  In  these  systems,  the  introduction  of  nitro  groups  also  has  an  overall  stabilizing  effect,  which  may  be 
due  to  further  delocalization  of  e-electrons,  reflecting  the  strong  electron-withdrawing  power  of  the  nitro  group. 


Introduction 

In  a  recent  computational  study,1  we  have  used  the  well-es¬ 
tablished  isodesmic  reaction  procedure2-3  to  show  that  the  presence 
of  the  nitrogens  in  the  diaza  acyclic  imidazoline  (I)  and  piperazine 
(II)  rings  confers  an  added  degree  of  stability  to  these  molecules. 


N^N 


I 


0 


'N' 

II 


We  attributed  this  to  ^-conjugation  of  the  nitrogen  lone  pairs.4-5 


(1)  Murray,  J.  S.;  Redfern,  P.  C.;  Lane,  P.;  Politzer.  P.;  Wilier.  R.  L.  J. 
Mol.  Struct.  iTHEOCHEM),  in  press 

(2)  Hehre,  W.  S.;  Ditchfield,  R.;  Radom,  L.;  Pople,  J.  A,  J.  Am.  Chem. 
Soc.  1970,  92,  4796. 

(3)  Hehre,  W.  J.;  Radom,  L.;  Schleyer,  P  v.  R.;  Pople,  J.  A.  Ab.  Initio 
Mnlerulnr  Orbital  Thcur),  Wiley.  New  York,  I7S6. 

(4)  Dewar,  M.  J.  S.  J.  Am.  Chem.  Soc.  1984,  !06,  669 

(5)  Minkin,  V.  I.;  Glukhovtsev,  M,  N.;  Simkin,  B.  Ya.  J.  Mol.  Struct. 
( THEOCHEM)  1988,  181.  93 


We  now  show  that  this  is  a  general  effect  that  applies  to  a  large 
group  of  aliphatic  and  aiicyclic  aza  systems  and  is  even  amplified 
by  the  presence  of  nitro  substituents.  These  stabilizing  tendencies 
are  of  particular  interest  because  of  the  importance  of  many  of 
these  molecules  as  high-energy  systems. 


Methods  and  Procedure 

The  isodesmic  reaction  procedure2-3  is  a  means  of  studying 
anomalous  energy  effects  in  molecules.  An  isodesmic  reaction 
is  a  hypothetical  chemical  process  in  which  the  number  of  bonds 
of  each  formal  type  remains  the  same  on  each  side  of  the  equation, 
but  their  mutual  relationships  are  changed.  Representative  iso¬ 
desmic  reactions  for  cyclopentane  (CSHI0)  and  imidazolidine 
(CjN2H8)  are  given  in  eqs  1  and  2,  respectively.  The  values 


5C2H6  —  C5H,0  +  5CH4  (1) 

4CH3NH2  +  C2H6  -  CjNiH,  +  3CH«  +  2NH,  (2) 


for  reactions  such  as  these  reveal  any  deviations  from  bond  energy 
additivity  and  arc  therefore  interpreted  as  being  due  to  special 
energy  effects  associated  with  the  molecule  being  investigated, 
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TABLE  I:  Comparison  of  Isodesmic  Reaction  Energies  Obtained  at 
Various  Computational  Levels 


kcal/mol 


molecule 

3-21G 

6-3 1 G* 6 

MP2/6-31G1 

CHjCH2CHj 

-1.4 

-0.8 

-2.2 

NH(CHj)2 

-2.3 

-2.2 

-4.1 

CHjCH2N02 

-7.2 

-3.5 

-5.0 

nh2ch2nh2 

-10.5 

-7.9 

-10.1 

N^N 

CJ 

-113 

-9.0 

-16.8 

*  Reaction  energies  are  obtained  from  energies  of  the  reactant  and 
product  molecules  optimized  at  the  respective  computational  levels. 
‘The  total  energies  of  these  molecules,  in  the  order  listed,  are 
-1 18.263  65,  -134.238  85.  -282.701  02,  -150.237  1  1,  and  -227.128  1 1 
hartrees.  cThe  total  energies  of  these  molecules,  in  the  order  listed,  are 
-1 18.67441,  -134.678  10,  -283.52021,  -150.70791,  and  -227.86269 
hartrees. 

e.g.,  strain,  resonance  stabilization,  steric  effects,  etc.  (A E  <  0 
indicates  the  presence  of  stabilizing  factors,  while  A£  >  0  implies 
destabilization.) 

We  have  used  the  ab  initio  self-consistent-field  (SCF)  molecular 
orbital  approach  (gaUSSIanss6)  to  optimize  molecular  geometries 
and  obtain  their  energies  at  the  3-2 1G  level.  SCF  methods  have 
been  shown  to  give  satisfactory  results  for  isodesmic  energy 
differences,3-1  when  viewed  on  a  relative  basis. 

As  part  of  our  analysis,  we  have  computed  the  electrostatic 
potentials  in  the  neighborhoods  of  the  aza  nitrogens  in  some  of 
the  molecules  studied.  The  electrostatic  potential  V(r)  that  is 
created  at  any  point  r  in  the  space  around  a  molecule  by  its  nuclei 
and  electrons  is  given  rigorously  by  eq  3.  ZA  is  the  charge  on 
nucleus  A,  located  at  RA,  and  p(r)  is  the  electronic  density  function 
of  the  molecule,  which  we  obtain  from  the  molecular  wave 
function. 


V(r) 


r-A 


d  r 


(3) 


V(r)  is  a  real  physical  property  that  can  be  determined  ex¬ 
perimentally  as  well  as  computationally8  and  is  well  established 
as  a  means  of  interpreting  and  predicting  chemical  reactive  be¬ 
havior,  especially  toward  electrophiles.8-10  In  the  present  work, 
however,  we  use  it  as  an  indication  of  the  degrees  of  delocalization 
of  nitrogen  lone  pairs.  Following  up  on  earlier  work,"  we  shall 
also  use  the  position  of  the  most  negative  value  of  the  potential 
as  a  means  for  defining  dihedral  angles  involving  the  nitrogen  lone 
pairs.  We  have  computed  V(r)  at  the  STO-5G  level,  using  the 
3-2 1G  optimized  structures.  Extensive  investigations  have  shown 
that  ab  initio  SCF  approaches  produce  generally  reliable  repre¬ 
sentations  of  K(r);9-,J'14  in  particular,  we  have  had  good  experience 
with  the  STO-5G  basis  set. 


Results 

In  order  to  investigate  the  dependence  of  isodesmic  reaction 
energies  upon  the  computational  level,  we  have  carried  out  a 


(6)  Frisch,  M.  J.;  Head-Gordon,  M.;  Schlegel,  H.  B.;  Raghavachari,  K,; 
Binkley,  J.  S.;  Gonzalez.  C.;  Defrees,  D.  J.;  Fox,  D.  J.;  Whiteside,  R.  A.; 
Sceger,  R.;  Melius,  C.  F.;  Baker,  J.;  Martin,  R.;  Kahn,  L.  R.;  Stewart,  J.  J. 
R;  Fluder,  E.  M.;  Topiol,  S.;  Pople,  J.  A.  GAUSSIAN  88,  Gaussian  Inc.: 
Pittsburgh,  PA,  1988. 

(7)  HinchlifTe,  A.  Computational  Quantum  Chemistry ;  Wiley:  New  York, 
1988;  Chapter  5, 

(8)  Politzer,  P.,  Truhlar,  D.  G.,  Eds.  Chemical  Applications  of  Atomic  and 
Molecular  Electrostatic  Potentials',  Plenum  Press:  New  York,  1981. 

(9)  Scrocco,  E.;  Tomasi,  J  Adv.  Quantum  Chem.  1978,  II,  115. 

(10)  Politzer,  P.;  Daiker,  K.  C.  In  The  Force  Concept  in  Chemistry:  Deb 
B.  M..  Ed.;  Van  Nostrand  '.cinhoiu.  York,  1981;  p  2V4, 

(11)  Murray,  ].  S.;  Sukumar,  N.;  Ranganathan,  S.;  Politzer,  P.  Ini.  J. 
Quantum  Chem.,  in  press. 

(12)  Scrocco,  E.;  Tomasi,  J.  Topics  in  Current  Chemistry,  No.  2; 
Springer- Verlag:  Berlin,  1973;  p  95. 

(13)  Ghio,  C.;  Tomasi,  J.  Theor.  Chim.  Acta  1973,  30,  151. 

(14)  Almloff,  J,;  Stogard,  A.  Chem.  Phys.  Lett.  1974,  29,  418. 


TABLE  II;  Calculated  3-21G  Molecular  Energies  and  Isodesmic 
Reaction  Energies 


molecule 

total  energy, 
hartrees 

isodesmic’ 

kcal/mol 

relative 

AE  isodesmic, 

kcal/mol 

ch4 

-39.976  88 

0.0 

c 

NH3 

-55.872  20 

0.0 

c 

c2h6 

-78.793  95 

0.0 

c 

CHjNH2 

-94.681  66 

0.0 

c 

c3h8 

-1 17.61330 

-1.4 

0.0 

NH(CHj)2 

-133.49485 

-2.3 

-0.9 

nh2ch2nh2 

-149.403  12 

-10.5 

-9.1 

NH2(CH2)2NH2 

-188.21757 

-8.8 

c 

NH2CH2NHCHj 

-188.21367 

-11.2 

c 

111* 

o 

-194.08055 

+3.0 

0.0 

IV 

-209.974  53 

-5.7 

-8.7 

1* 

N^N 

-225.86361 

-11.3 

-14.3 

w 

V* 

O 

-232.91691 

-9.1 

0.0 

VI 

o 

-248.803  26 

-13.0 

-3.9 

11* 

CN) 

-264.689  25 

-16  7 

-7.6 

Vli 

-264.694  57 

-20.1 

-11.0 

G 

VIII 

-303.505  64 

-16  3 

c 

G 

IX 

-372.962  58 

-41.1 

c 

CG 

X 

-450.602  21 

-44.5 

c 

CO 


“Relative  isodesmic  energies  are  referred  to  the  corresponding  hy¬ 
drocarbon.  ‘  Data  are  taken  from  ret  I.  ‘  Not  applicable. 


TABLE  III:  Calculated  3-21G  Molecular  Energies  and  Isodesmic 
Reaction  Energies  for  Some  Nitro  Derivatives 


total  energy.  A£ilotem,c, 
molecule  hartrees  kcal/mol 


ch,no2 

-242.255  85 

0.0 

nh2no2 

-258.13800 

0.0 

ch3ch2no2 

-281.084  40 

-7.2 

(CHj)2NN02 

-335.77489 

-11.3 

N02NH(CH2)2NHN02 

-592.75641 

-13.1 

XI 

0,N  NO, 

N  ^  ✓ 

N  N 

\_y 

-630.413  58 

-22.9 

XII 

no2 

-669.244  36 

-31.4 

I 

NOj 


comparison  of  3-21G,  6-31G*,  and  MP2/6-31G*  results  for  a 
limited  number  of  systems.  MP2/6-31G*  geometry  optimizations 
for  nitroethane  and  imidazolidinc  were  carried  out  using  the 
Cambridge  Analytic  Derivatives  Package  (cadpac)15  at  the  San 
Diego  Supercomputer  Center;  GAUSSIANSS*  was  used  for  all  other 
computations.  Our  results  are  given  in  Table  I.  For  each  of  these 
three  computational  approaches,  the  A E  values  show  the  same 
trend,  becoming  more  negative  in  going  from  propane  to  imida- 


(15)  Amos,  R.  D.;  Rice,  J  E.;  Handy,  N  C.  The  Cambridge  Analytic 
Derivatives  Package,  1987. 
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zolidine.  These  data  indicate  that  isodesmic  reaction  energies 
computed  at  different  levels  will  lead  to  similar  conclusions  re¬ 
garding  anomalous  energy  effects  when  viewed  on  a  relative  basis. 
For  the  present  study  we  have  chosen  the  3-2 1G  approach  as  being 
computationally  the  most  practical. 

Our  calculated  3-2 1G  molecular  and  isodesmic  reaction  energies 
for  a  group  of  aliphatic  and  alicyclic  systems  and  some  of  their 
nitro  derivatives  are  presented  in  Tables  II  and  III,  respectively. 
Positive  A£  values  are  commonly  interpreted  as  indicating  strain. 
For  example,  the  3-2 1G  isodesmic  reaction  energy  of  cyclopropane 
is  31  kcal/mol,16  which  is  quite  close  to  its  widely  quoted  strain 
energy  of  27  kcal/mol.17  The  AE  value  of  cyclopentane,  +3 
kcal/mol  (Table  II),  indicates  a  low  degree  of  strain,  in  agreement 
with  its  reported  strain  energy  of  7  kcal/mol.17  All  other  molecules 
in  Table  II  have  negative  isodesmic  energies,  reflecting  varying 
degrees  of  stabilization,  which  we  attribute  to  cr-conjugation.4-5 

Looking  first  at  Table  II,  it  is  important  to  note  that  the  sat¬ 
urated  unstrained  hydrocarbons  propane  and  cyclohexane  have 
small  negative  calculated  AE  values.  These  results  are  consistent 
with  theoretical  studies  that  have  shown  that  a  stabilizing  delo¬ 
calization  of  (r-electrons  occurs  in  saturated  hydrocarbons.1819  The 
aza  systems  have  more  negative  AE'%  than  their  hydrocarbon 
analogues,  presumably  reflecting  the  additional  presence  of  ni¬ 
trogen  lone  pairs  that  are  delocalized  to  some  degree  through  the 
<7-bond  framework.  This  extra  stabilization  occurs  to  a  greater 
extent  in  the  cyclic  than  in  the  noncyclic  systems  and  increases 
even  more  strikingly  when  there  are  fused  rings.  Within  any 
particular  aliphatic  or  alicyclic  series  of  molecules  with  the  same 
number  of  first-row  atoms,  the  stabilization  is  greater  as  the 
number  of  nitrogens  increases.  This  point  is  brought  out  even  more 
strikingly  by  the  relative  A£isod<sm,c  values  that  are  given  in  Table 
II. 

It  is  interesting  to  observe  that  the  1,3-diaza  systems  show  a 
greater  degree  of  stabilization  than  do  their  1,4-diaza  isomers. 
This  can  be  seen  in  the  AE  values  for  the  pair  NH2CH2NHCH3 
and  NH2(CH2)2NH2  and  again  for  1,3-  and  1,4-diazacyclohexane 
(VII  and  II)  [Table  II].  The  fact  that  the  fused-ring  systems  IX 
and  X  have  large  calculated  AE  values  (significantly  more  than 
twice  the  AE's  of  I  and  II)  may  be  due  to  the  additional  two 
1 ,3-diaza  groupings  that  occur  in  each  by  virtue  of  being  fused. 

Isodesmic  reaction  energies  for  nitro  derivatives  of  ethane, 
dimethylamine,  1,2-diaminoethanc,  imidazoline,  and  piperazine 
are  listed  in  Table  III.  The  AE  values  for  these  nitro  systems  are 
larger  in  magnitude  than  those  for  the  corresponding  unsubstituted 
molecules,  indicating  that  the  introduction  of  nitro  groups  has  an 
overall  stabilizing  effect  upon  these  systems.  This  additional 
stabilization  may  be  due  to  enhanced  delocalization  of  the  a- 
electrons,  particularly  the  amine  lone  pairs,  resulting  from  the 
strong  electron-withdrawing  power  of  the  nitro  group.  (Nitro 
groups  do  not  invariably  have  a  stabilizing  effect;  when  substituted 
on  tertiary  carbons  in  strained  molecules,  they  are  found  to  de¬ 
stabilize  the  systems.20) 

Discussion 

We  have  suggested  that  the  extra  stability  provided  by  the 
replacement  of  carbons  by  nitrogens  in  aza  aliphatic  and  alicyclic 
systems,  sucii  as  those  in  Table  II,  is  largely  due  to  lone-pair 
cr-conjugation.  As  a  measure  of  the  extent  of  this  conjugation, 
we  examine  the  electrostatic  potentials  V(r)  in  the  vicinities  of 
the  aza  nitrogens. 

C(r)  for  aliphatic  and  aromatic  amines,1-21-26  hydrazines,11  and 


(16)  Politzer,  P.;  Lane,  P.  Struct.  Chem.,  in  press. 

(17)  Greenberg  A.;  Liebman,  J,  F.  Strained  Organ, c  Minecules.  Auuicrr.,. 
Press:  New  York,  1978 

(18)  Herndon,  W,  C.  J.  J.  Chem.  Educ.  1979,  56,  448. 

(19)  Herndon,  W.  C.  J.  Tetrahedron  Lett.  1979,  3801. 

(20)  Murray,  J.  S.;  Seminario,  J.  M.;  Politzer,  P.  Unpublished  work 

(21)  Politzer,  P.;  Abrahmsen,  L ;  Sjoberg,  P.  J.  Am.  Chem.  Soc.  1984,  106, 
855. 

(22)  Politzer,  P.;  Bar-Adon,  R,  J.  Phys.  Chem.  1987,  91.  3069. 

(23)  Politzer,  P  ;  Elminyawi,  I.  M.;  Lane,  P,;  McKenney,  R.  L.,  Jr.  J.  Mol. 
Struct.  ( THEOCHEM ),  in  press 

(24)  Murray,  J.  S.;  Politzer,  P  Chem.  Phys.  Lett.  1988,  152.  364. 


Figure  1.  Calculated  electrostatic  potential  of  1,3-diazacyclohexane 
(  VII),  in  kilocalories  per  mole,  in  a  plane  passing  through  the  two  NH 
groups.  The  projections  of  other  nuclear  positions  are  shown  by  their 
atomic  symbols.  Dashed  contours  correspond  to  negative  potentials.  The 
positions  of  the  most  negative  potentials  are  indicated;  the  value  is  (•) 
-72.0. 

TABLE  IV:  Calculated  Electrostatic  Potential  Minima  ( V'^N)]  of 
Aza  Nitrogens 


FmJN). 

molecule  kcal/mol 


NH3 

NH(CH3)2 

-109.3“ 

-100.0 

IV 

N 

o 

-97.2 

VI 

0 

-95.5 

1 

N^N 

-95.3'’ 
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3  is  taken  from  ref  24. 

6 The  y„,„( N )  data  for 

imida70lidine  and  piperazine  are  taken  from  ref  I. 

nitrogcn-conlaining  heterocycles l0-27-28  are  largely  dominated  by 
strong  negative  regions  attributed  to  the  nitrogen  lone  pairs, 
indicative  of  sites  attractive  to  electrophiles.  An  evnrmle  of  these 
characteristic  poiemials  is  given  by  the  k(e)  of  1 ,3-diazacyclo- 
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hexane  (VII),  shown  in  Figure  I ;  each  ring  nitrogen  has  a  negative 
potential  region  with  a  most  negative  value,  Fmin(N),  associated 
with  it.  The  Fmjn(N)  values  for  VII  and  various  other  aza  systems 
are  listed  in  Table  (V.  [Fmin(N)  for  VII  is  not  located  in  the  plane 
shown  in  Figure  1.] 

In  an  earlier  study,  we  have  shown  that  the  Fmin(N)  values  of 
NH2-X  molecules  reflect  the  total  electron-attracting  tendencies 
of  the  substituents  X  24  The  decrease  in  magnitude  of  Fmj„(N) 
relative  to  NH3  is  viewed  as  an  indication  of  the  degree  to  which 
the  nitrogen  lone  pair  electronic  density  is  redistributed  or 
“delocalized”  throughout  the  remainder  of  the  molecule.  We 
interpret  the  Kmin(N)  values  in  Table  IV  in  a  similar  manner. 

In  general,  these  Fmi„(N)  decrease  in  magnitude  as  the  iso- 
desmic  reaction  energies  of  the  corresponding  molecules  become 
more  negative.  (This  trend  holds  rigorously  within  the  five-  and 
six-membered  cyclic  systems  taken  separately.)  Thus,  increased 
delocalization  of  the  nitrogen  lone  pairs,  as  measured  by  decreases 
in  the  magnitudes  of  the  Kmjn(N)  values,  roughly  parallels  an 
increase  in  stability.  These  data  support  the  view  that  the  en¬ 
hanced  stability  in  aliphatic  and  alicyclic  aza  systems  is  linked 
to  lone-pair  ^-conjugation. 

It  has  been  pointed  out  that  the  1,3-diaza  molecules,  NH2C- 
H2NHCH3  and  1,3-diazacyclohexane,  have  a  greater  degree  of 
stabilization  than  their  1,4-diaza  isomers  (Table  II).  We  suggest 
that  this  extra  stabilization  comes  about  because  these  systems 
can  assume  conformations  that  permit  the  <r-delocalization  of  each 
nitrogen  lone  pair  to  reinforce  that  of  the  other.  The  tendency 
to  assume  such  geometries  has  been  noted  and  studied  earlier4-29"31 
and  has  been  called  the  anomeric  effect.  It  is  found  in  R'-X- 
C-Y-R  systems,  where  X  and  Y  are  heteroatoms  with  lone  pairs, 
that  the  dihedral  angle  defined  by  the  lone  pair  of  X  and  the  X, 
C,  and  Y  nuclei  tends  to  be  close  to  130°.  This  has  been  ra¬ 
tionalized  in  the  molecular  orbital  framework  as  stabilization 
occurring  via  delocalization  of  electronic  density  from  a  non¬ 
bonding  molecular  orbital  of  X  into  the  <r*-orbital  of  the  adjacent 
C-Y  bond.30-31 

Recent  ab  initio  3-21G  calculations  of  NH2CFI2NH2,  the 
smallest  molecule  containing  the  1,3-diaza  moiety,  have  reported 
that  its  most  stable  conformer  possesses  a  calculated  structure 
consistent  with  the  anomeric  effect;32  the  lone-pair-N-C-N  di¬ 
hedral  angles  are  180°.  (We  have  included  NH2CH2NH2  in 
Tables  I  and  II.)  The  1,3-diaza  rings  I,  VII,  and  VIII  contain 
this  methylenediamine  unit  and  have  calculated  F^NJ-N-C-N 
dihedral  angles  that  increase  from  129°  in  I  to  169°  in  VII  to 
174°  in  VIII.  This  is  paralleled  by  a  decrease  in  the  N-C  bond 
distances  within  the  N-C-N  groupings  of  I,  VII,  and  VIII;  these 
are  1.486,  1468,  and  1.459  A,  respectively.  These  data  are 
consistent  with  the  delocalization  of  the  nitrogen  lone  pairs  in¬ 
creasing  as  the  structural  criteria  for  the  anomeric  effect  are  more 
fully  satisfied.  It  is  interesting  to  note  that  the  Kmi„(N)  values 
decrease  in  magnitude  in  the  same  direction,  from  the  five-  to  the 
seven-membered  1,3-diaza  rings,  which  is  also  consistent  with 
increased  nitrogen  lone  pair  cr-conjugat-on.  That  the  Fmin(N)- 
N-C-N  dihedral  angles  in  I  deviate  significantly  from  180°  may 
be  a  result  of  a  lack  of  flexibility  in  the  imidazolidine  ring,  which 
prevents  a  more  preferred  conformation  with  regard  to  the  position 
of  the  nitrogen  lone  pairs  and  the  adjacent  C-N  bonds. 

We  find  that  the  most  stable  conformers  of  the  fused  ring 
systems,  IX  and  X,  have  geometries  that  preclude  the  anomeric 
effect  from  operating  in  both  directions  within  their  respective 
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I, 3-diaza  groupings.  For  example,  our  calculated  Fmjn(N')- 
N'-C-N  and  F?m(N)-N-C-N'  dihedral  angles  in  X  are  176° 
and  61°,  respectively.  Furthermore,  the  Fmin(N')  and  the  N'-C 
bond  distances  (-88.9  keal/mol  and  1.455  A)  are  smaller  in 
magnitude  than  the  Fmjn(N)  and  the  N-C  bond  distances  (-98.5 
keal/mol  and  1 .469  A).  These  data  are  consistent  with  lone-pair 
^-conjugation  occurring  to  a  greater  extent  for  the  nitrogens, 
designated  as  N'  in  Table  IV,  which  can  assume  conformations 
that  closely  satisfy  the  requirements  for  the  anomeric  effect.  A 
similar  situation  is  observed  for  the  fused  five-membered  ring 
system,  IX,  which  has  four  1,3-diaza  groupings  in  succession 
around  the  outer  framework  of  the  molecule.  In  one  direction 
around  the  ring,  the  Kmjn(N)-N-C-N  dihedral  angles  are  150° 
or  151°;  in  the  opposite  direction,  these  angles  are  alternately  83° 
and  26°.  That  each  nitrogen  in  IX  satisfies  the  anomeric  effect 
criterion  in  one  direction  is  reflected  in  the  similar  Kmi„(N)  and 
Fmin(N')  values  (Table  IV). 

Stabilization  due  to  ring  nitrogens  in  alicyclic  systems  may  be 
a  key  factor  involved  in  understanding  some  of  the  anomalies  in 
the  thermolysis  patterns  of  nitramines.  It  has  been  found  in  the 
fast  thermolysis  of  many  secondary  nitramines  that  the  relative 
initial  concentration  of  N02  depends  qualitatively  on  the  N-N 
bond  distance  and  the  average  asymmetric  stretching  frequency 
of-N02  MN02)].33  Molecules  with  long  N-N  bonds  and  high 
vas(N02)  values  are  often  strong  liberators  of  N02  upon  fast 
thermolysis;  those  with  short  N-N  bonds  and  low  stretching 
frequencies  often  do  not  liberate  N02  or  do  so  in  competition  with 
other  reactions.  1 ,3-Dinitroimidazolidine  (XI)  deviates  from  this 
qualitative  trend;  it  has  a  relatively  short  calculated  N-N  bond 
length  ( 1 .344  A)1  but  liberates  N02  by  cleavage  of  the  N-N  bond 
during  thermolysis.33  This  exceptional  behavior  is  believed  to 
reflect  an  anomalously  high  stability  of  the  imidazolidine  ring,34 
consistent  with  the  stabilization  shown  for  this  molecule  in  Table 

II. 

Summary 

The  isodesmic  reaction  procedure  has  been  used  to  study  the 
effect  of  replacing  carbons  by  nitrogens  in  a  group  of  aliphatic 
and  alicyclic  systems.  In  all  instances,  we  find  the  aza  systems 
to  possess  significant  degrees  of  added  stability,  which  we  attribute 
to  cr-conjugation.  This  presumably  involves  the  nitrogen  lone  pairs, 
which  are  delocalized  to  some  extent  through  the  <r-bond 
framework.  The  stabilization  due  to  the  presence  of  nitrogens 
is  greater  in  cyclic  systems,  particularly  when  rings  are  fused,  and 
it  increases  with  the  number  of  nitrogens.  Decreases  in  the 
magnitude  of  the  most  negative  electrostatic  potentials  of  the 
nitrogens,  which  can  be  viewed  as  a  measure  of  increasing  delo¬ 
calization  of  the  nitrogen  lone  pairs,  roughly  parallel  the  increase 
in  stability.  Isodesmic  reaction  energies  for  nitro  derivatives 
indicate  that  in  these  systems  the  introduction  of  nitro  groups  also 
has  an  overall  stabilizing  effect,  which  may  be  a  consequence  of 
further  delocalization  of  ^-electrons,  particularly  the  amine  lone 
pairs,  due  to  the  strong  electron-withdrawing  effect  of  the  nitro 
group. 
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